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We study the linearized equations describing the propagation of gravitational waves 
through dust. In the leading order of the WKB approximation, dust behaves as a non- 
dispersive, non-dissipative medium. Taking advantage of these features, we explore the 
possibility that a gravitational wave from a distant source gets trapped by the gravitational 
field of a long filament of galaxies of the kind seen in the large scale structure of the 
Universe. Such a waveguiding effect may lead to a huge magnification of the radiation 
fiux from distant sources, thus lowering the sensitivity threshold required for a successful 
detection of gravitational waves by detectors like VIRGO, LIGO and LISA. 



Abstract 



PACS: 04.25, 04.30-w, 98.62.Sb 



1 Introduction 



Detecting gravitational waves is one of the greatest challenges for modern physics both from 
the theoretical and the experimental points of view ||, ^. The interest of this enterprise 
arises from the fact that, in case of success, it will open an entirely new window on the Universe 
and, as our past experience with radio astronomy teaches, this might have a profound impact 
on our vision of the Universe. Unfortunatly, the detection of a gravitational wave has turned 
out to be a very difficult task and as of now, almost fourty years after the early pioneering work 
of Weber no one has succeded yet. The main source of the difficulties is the fact that most 
of the known astrophysical sources should radiate exceedingly small amounts of gravitational 
radiation. The strongest sources should be massive, highly non-sperical systems with large 
internal kinetic energies like colliding and coalescing black holes and neutron stars. Sources 
of this type last for a short time and should also be very rare, with the consequence that it will 
be necessary to search large regions of the Universe in order to have a reasonable event rate. 
For sources in the high frequency band (1-^ 10^ Hz), accessible to ground based interferometers 
like VIRGO and LIGO, a rate of several coalescence events per year requires extending the 
search at distances from the Earth of the order of 10^ 10^ Mpc. At such large distances 
even the dimensionless amplitudes h of the waves from such strong sources should be in the 
range of /i ~ 10~^^ -i- 10~^'^, which still represents an extremely small signal. There are hopes 
that the LIGO/VIRGO collaborations |||],||2| will be able to reach or at least closely approach 
the sensitivity needed for the detection of these signals. The same considerations apply to 
the sources in the low- frequency band (10~^ -i- 1 Hz) accessible to the forthcoming orbiting 
LISA detector. Examples of such sources from distant galaxies, at distances of the order of 
one Gpc from the Earth, are coalescing massive black holes, or compact bodies spiraling into 
one of them. For a more thorough discussion of these issues and an updated guide to the 
bibliography we refer the reader to the recent paper by Thorne ^J. 

This being the state of things, it is our opinion that any mechanism leading to an am- 
plification of the signal from a cosmological source of gravitational waves, thus making its 
detection easier, should be very welcome. An example of such a mechanism, well established 
by now in the realm of optical and radio-wave astronomy, is provided by gravitational lens- 
ing. In fact, the possibility that the gravitational radiation from distant binary neutron stars 
may get magnified due to microlensing by compact bodies has been recently considered in 
the literature ||7|. In this paper we propose a different lensing eff'ect: we shall argue that the 
large scale galactic gravitational fields generated by the filamentary distributions of galaxies 
seen in the large scale structure of the Universe, primordial cosmic strings evolved in today 
observable structures, may act as waveguides for gravitational waves. Since a gravitational 
wave interacts very weakly with practically all kinds of matter, a wave trapped in the interior 
of such a filament might be able to traverse cosmological distances without being absorbed or 
scattered in a significant way, and most importantly without undergoing the 1/R attenuation 
that charaterizes propagation in vacuum. If a significant fraction of the radiation from distant 
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sources is waveguided to us in this way, this might lower in an appreciable way the minimum 
sensitivity needed for detecting the signal. 

The idea of gravitational conductors, waveguides and circuits is not new per se. For ex- 
ample, it was considered in |^] where it was based on the hypothetical existence of a material 
having a large shear (the "respondium") and capable of reflecting a gravitational wave. Our ap- 
proach is completely different and not so exotic. Starting from standard General Relativity, we 
show that a filamentary distribution of dust is capable of trapping in its interior a gravitational 
wave: this result will be derived below by applying the (lowest-order) WKB approximation to 
the linearized Einstein's equations describing the propagation of a weak gravitational wave in 
a region of space filled with dust. Within these approximations, the trapping of a gravitational 
wave appears completely analogous to that of a ray of light inside a guide with a specific re- 
fraction index profile. Being based on the ray approximation, our approach cannot take into 
account interference effects, which should be negligible in the high frequency band accessible 
to VIRGO/LIGO as well as in the low frequency band accessible to LISA. In case they could 
not be neglected, one could resort to an approximation of the type of the Fock-Leontovich 
paraxial approximation Q, which, while using the short-wave approximation, allows also to 
describe interference. 

A waveguiding effect by cosmological structures analogous to that described in this paper 



was shown to be possible in Ref. pQ] for the case of electromagnetic waves and in Refs.|0, 12] 



for massless particles, using the Fock-Leontovich approximation. In Ref.|lC], in particular, this 
mechanism was used to explain in a simple way the huge luminosity of quasars, as compared 
to their high redshift. They would not be exotic objects hiding in their interior some unknown 
engine that releases enormous amounts of energy; rather, they might just be protogalaxies 
whose luminosity is preserved by a waveguide. Furthermore, certain gravitational lensing 



effects such as "twin" and "brother" quasars [13| could be explained in a similar manner. 



The paper is organized as follows. In Sect. 2 we present the linearized theory of grav- 
itational radiation in a region of space-time filled with dust. We show that, in the leading 
WKB approximation, dust behaves as a non-dispersive and non-dissipative medium. In Sect. 
3, we use this approximation to study the possibility that a gravitational wave, coming from 
a distant source, could be trapped by the gravitational field of a long filamentary structure, 
thus yielding to a waveguide effect. We estimate the focalization length, the amount of lumi- 
nosity preserved, and the effective feasibility of the effect. A discussion of the results and the 
conclusions are drawn in Sect. 4. 
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2 Linearized Theory 



In this Section, we shall derive the linearized equations describing the propagation of a weak 
gravitational wave through a region of space-time filled with dust. 0. The same problem was 
discussed in Ref.|16], but the equations derived below are simpler than those of p^ . Before 
the passage of the wave we thus have an unperturbed background {OabTUa, p), where Qab is the 
space-time metric while Ua and p denote, respectively, the field of velocities and the density of 
the dust. The background is assumed to satisfy the exact Einstein's equations in the presence 
of dust: 

Gab = KpUaUb , (1) 

where k = SvrG/c^, G being Newton's constant, while Gab is the Einstein's tensor 

Gab = Rab — -^9abR ■ (2) 

We also have 

p > , (3) 

while 

UaU'' = -1 . (4) 

We recall that the conservation law for the stress-energy of the dust, combined with the 
normalization condition (^) implies that the velocities u"" satisfy the geodesic equation: 

u^V^u" = , (5) 

while the density p satisfies the conservation equation: 

^aipu'') = . (6) 

We now imagine perturbing the background solution by a small perturbation {gab + {jab^^a + 
Ua, p+p)- Plugging the perturbed solution into Einstein's equations and retaining all the terms 
up to first order in the perturbation, we get the following linearized equations: 



where 



Gab = K[pUaUb + piUaUb + UaUb)] , (7) 



Gab= - ^V^Vc7afe + V(„V^7b)c - ^SafeV^V^Tcd 



1^- 1 ^cd. 



2^ 

+ R'ab "led + R(^a%c - \Rlab + ^QabR'^'lcd ■ (8) 

In this equation Va and Rahc represent, respectively, the covariant derivative and the Riemann 
tensor relative to the background metric gab- Indices are raised and lowered using gab and its 
inverse, while 7^6 is defined as 

(jab = Jab - -9ab 7 , (9) 



*We shall follow in this paper the notations and conventions of |l5|. In particular, the signature of the metric 
is chosen to be ( — h ++)• Moreover yl(bc) := ^{Atc + ^cb)- 
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where 



1 = 9 lab 



(10) 



We have to add to Eq. (^) the equation obtained hnearizing the analogue of the normahzation 
condition Eq. (^) for the perturbed velocity field Ua + Ua'- 



In deriving this equation, we used: 



and 



ac bd': 

-g 9 9cd 



u''=9'''ub-9'''n'9bc 



(11) 



fl2) 



(13) 



Thus, Eqs. (|^) and (11) are the equations that describe the propagation of a weak gravitational 
wave in a medium made of dust. We refer the reader to Ref. [16| for a proof of their linear 
stability. Let now be the projector onto the hyperplane orthogonal to the 4- velocity u": 



(14) 



Upon multiplying Eq. (0) by /le^j, we get: 

hlh) Gafe = . (15) 
Upon multiplying now Eq. (^) by and using Eq. ( p]l| ) to eliminate u°"Ua, we also get: 

Kp = Gabu"'v!' + KpgabU^'u^ . (16) 

Finally, by multiplying Eq. (|^ by and using again Eq. (0) together with Eq. (jl|), we get: 

KpUa = -{Gcdu'^u'^ Ua + GabU^) " '^P (Jcdu'^u'^ Ua ■ (17) 

Equations (|l5|), (p^) and (|l^ are completely equivalent to Eqs. (|^) and (11). This can be seen 
by plugging Eq. ( [l^ ) in Eq. ( p^ ) and then taking the terms involving the velocities Ua on the 
right hand side: upon using Eqs. (^) and ([T^), one then recovers the right hand side of Eq. 
(0). Finally, Eq. (11) can be obtained back from Eq. (|T7| ) by simply multiplying it by u"". 

Now, Eq. (|l^) is the only evolution equation governing the perturbation, while Eqs. (16) 
and (17) just represent constraints. It is quite remarkable that the evolution equations depend 
only on the perturbation of the metric, while the constraints are explicitly solved with respect 
to the perturbations of the velocities and density. Thus, we focus on Eqs. (15): they can be 
simplified a bit by a suitable choice of gauge. Under an infinitesimal diffeomorphism 



'Jab transforms as: 



lab = lab + 2V(„^b) - gab^^Cc ■ 



(18) 
(19) 
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Using the method discussed in 1 16 1, it is easy to prove that by means of a suitable choice of 
it is always possible to achieve the so-called transverse gauge: 



y'lab = . (20) 

We notice that this transversality condition does not fix completely the gauge, as it is preserved 
by any infinitesimal diffeomorphism such that: 

VVaCb + Rb'^a = . (21) 

In the transverse gauge, the second and third terms in the expression of Gab in Eq. (^) drop 
out and thus we are left with the following equations: 

/I 11 

'^Clhb / V7CV7 ^, . DC (i- . D C- , I . ^ . ncd' 



Kh) [^f^'-^clab - R'ab "led - R^:ib)c + l^Rlab " ^QabR'^ledj = , (22) 

that are to be solved together with Eq. (pO|). 

A further simplification arises if the wavelength A of the perturbation is much smaller than 
all the typical lengths that characterize the background. We thus assume that the distances 
over which the velocities Ua vary in an appreciable way and the average radius of curvature 
of the background \Rabcd\~'^^'^ ars both much larger than A. Under these conditions, it is 
legitimate to make the WKB ansatz: 

lab{x,e) = Re{/ai,(x,e) exp(i/e S{x))} , (23) 

where fab{x, e) is a slowly varying amplitude having the asymptotic expansion (for e — > 0): 



A,(x,e) = ^(^) /^:^(x). (24) 

n=0 



Plugging this expansion in Eq. ( p^ ) and equating to zero the coefficients of each power of e. 

An) 
lab 



one gets a system of recursive equations for the phase S{x) and for the amplitudes fab\x). To 



order e we get: 

(<7"'U.)(/i^/i^/i?) = 0, (25) 

where la is the wave vector: 

la := daS . (26) 
As for the transversality condition Eq. (^), it gives to order e~^: 

/Vi? = . (27) 

Equation (|2^) admits two types of solutions, depending on which of the two factors between 
the brackets vanishes. Consider first: 

Khif^^^ = . (28) 
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This equation implies that is of the form: 

= '^(a^b) ■= ]j{VaUb + UaVb) , (29) 

where Va in an arbitrary vector field. Upon multiplying both sides of the above equation 
by 2/^ and imposing Eq. (p7|), one finds that Va must be proportional to Ua, Va = vua- 
But then, multiplying Eq. (pg|) by one gets that fj^^ can be different from zero only if 
CO := —I'^Ua = 0. This means that the wavefronts are time-like surfaces, and then this type of 
solutions cannot describe a gravitational wave. For this reason, we shall not consider anymore 
this case in the rest of the paper. 



The second class of solutions of Eq. ( pS]) is characterized by phase-functions S{x) satisfying 
the eikonal equation: 

g'^HdaS) {dbS) = . (30) 

The wavefronts are thus null hypersurfaces (we exclude the possibility of singular wavefronts 
and thus assume that la does not vanish at any point) describing waves that propagate through 
the dust with the speed of light. We recall that the integral curves x"'{v), called rays, of the 
null vector field la- 

— = {g-'dtS){x{v)), (31) 

with S{x) a solution of the eikonal equation, are null geodesies of the background metric gab 
providing bicharacteristics of Eq. ([2^). 

By looking at the next order in the formal expansion of Eq. (22), at order one gets a 
continuity equation for the amplitude fj^l\x): 

KhtiVi + e)f^S = , (32) 



where we have set 



Now, Eq. (|32D is equivalent to 



:= rVa , e := ^Val" . (33) 



(V; + e)/i?=^;(a«b) , (34) 



where Va{x) is again an arbitrary vector field. Eq. (27) together with the geodesic condition 
satisfied by la, Vila = 0, imply that Va must be identically zero. To see this, multiply both 
sides of Eq. ( p4| ) by 2/". The left-hand side is then seen to vanish: 

2/'(V/ + 0)/l? = 2VK/Vi?) = O, 

and thus we have for the right-hand side: 

= 2l^V(^aUb) = {l^Vb)Ua - UJ Va . (35) 

This equation implies that Va must be proportional to «„: Va = v Ua', but then it easily follows 
that it must vanish because after multiplying Eq. (^) by l°'l'^ we get 
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which imphes that v vanishes, because w, being the time-component of a non-vanishing null 
vector, is different from zero. We conclude that the amplitude fj^^ satisfies the transport 
equation: 

(V« + ^)/i? = 0. (36) 

It follows from this equation that the expansion of Gab vanishes up to order included and 
thus, according to Eqs. (16) and ([TTD, the dust suffers no density or velocity perturbations to 
this order. 



As it is usual in the WKB approximation, Eq. (|36| ) converts into an ordinary first-order 
differential equation for the restriction fj^\x{v)) of the amplitude fj^\x) to the rays ( |3l| ) and 



thus its solutions are uniquely determined by the knowledge of on a surface S cutting each 
ray at one point. We can further simplify Eq. p^) by a suitable choice of a tetrad field, which 



we now describe. We take and k°' := h'^l^ as the first two elements of the tetrad, while the 
remaining two, that we call e^, are defined as follows. We pick up an arbitrary point on 
each ray and choose for e^, there two arbitrary space-like unit vectors such that: 

eX = eX = 0, e'-el = 0. (37) 

The vectors then transported along the ray according to the "quasiparallel" transport 

rule [ p!7| ] given by the equation 

V;e*« = -e'^Viqt , (38) 
where is the orthogonal projector onto the plane spanned by u"" and A;": 

qt = -u^Ub + ^k^kb . (39) 

In [0] it is shown that this rule of transport preserves inner products and that remain 
orthogonal to u"" and k"" all along the ray. Using the tetrad so defined, one can see that the 
most general amplitude fj^^ satisfying the gauge condition Eq. (p7|) can be decomposed as 

/i? =A+e+, + A^ el, + Atre'^, + iJ.ej J^) + Blah , (40) 



where 



a„i „2„2 X _plp2 , 2 1 

^ab ~ *^a*^b ^ ^a^b ' 

tr I A , „2„2 



C = e^ei + eiet . (41) 

e^j, and e^j, represent the usual traceless transverse (T-T) modes of a gravitational wave, while 
e**^ is a transverse mode with non- vanishing trace; as for the last two terms in Eq. (40), 



they represent longitudinal modes. The coefficients A^, Ax, An-, Bi and B satisfy a simple 
transport equation. Upon inserting Eq. (Hol) in Eq. (B^) and taking its trace, we find: 



(V; + e)Atr = . (42) 
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Contracting now Eq. (36) with any of the traceless modes appearing in Eq. (|40|), in view of 
the fact that they are orthogonal to each other, and observing that e^V/e-'" = 0, which is a 
consequence of the transport rule Eq. (38), one obtains: 



{Vi + e)B, = {Vi + 9)B =0 



(43) 



We now show that by means of a gauge transformation of the type (21) it is always possible to 
achieve that Atr, Bi, and B vanish on some initial surface S and, due to the above transport 
equations, this implies that they vanish all along the rays crossing S. For this purpose, we 
consider a family of vector fields ^(x, e) of the form 



n=l ^ 



(X) . 



(44) 



This ensures that gauge transforming a jab of the form ( 23 ) according to Eq. ( |l9|) the new jab 

AO) 

lab 

P(0)/ _ ^(0) f(l) _ ^,c.(l)^ 



is still of the form (23). In particular, we find for fjf'^ 



a' = a+^hae,)'-m^)9ab 



(45) 



Upon imposing Eq. ( pl| ) on the field ^(x, e), one can see that the coefficient ^i^^(x) must satisfy 
the transport equation 

(Vi + 0)d')=O. (46) 

The important thing is that the value of ^i^^ on S is completely arbitrary and it is easy to 
verify that it always possible to choose it in such a way that 

= B[ = B' = {) on S . (47) 



Consider now the effective energy-momentum tensor T"^ [18| and the effective "graviton num- 
ber" A^"^, respectively, defined as 



T 

N' 



ah 



l^+l^ + l^xP) I'^l'' , 



4k 
1 



(|A+|^ + |^xP)r . 



(48) 



The transport equations for A^ and together with the geodesic equation for la, imply that 
their covariant divergences vanish: 



Vaf"^ = 



VaN" = . 



(49) 



Using the identity F^^ = dag/ (2^), where g =det || gab || , we can rewrite the equation VaN"- = 
in the form of a conservation law: 



1 



4KWaiV° = dai\A\' n + \A\'r-dag 

29 



1 



2ia\ 



-.da{^\A\'l 







(50) 
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where A'^ = + |AxP)- Upon integrating this equation on the 4-volume spanned by a 

ray bundle, and using Gauss' theorem, one proves that the number of gravitons in the bundle 
is conserved and observer-independent. 

The conclusion of this analysis is that, in the lowest WKB approximation, the physical 
degrees of freedom of a gravitational wave propagating through dust are represented, as in 
vacuum, by two traceless-transverse modes, e^^, and e^^. The transport equations for 
and Ax, Eqs. (^3|), show that these modes propagate along null geodesies of the background 
gravitational field, again as it happens in vacuum. The dust behaves as a non-dispersive, non- 
dissipative medium that preserves the polarization of the wave and the total graviton number 
in a ray bundle. 

3 The Waveguiding Effect 

It is well known that, when the WKB approximation is valid, a weak gravitational wave 
propagates in vacuum much like an electromagnetic wave. One is thus led to the expectation 
that a gravitational wave passing by some massive body may undergo a lensing effect, much 
in the same way as it happens to light. The results of the previous Section show that under 
many respects a gravitational wave propagating through dust does not behave differently than 
in vacuum. In particular, we saw that dust behaves as a non-dissipative medium which implies 
that a gravitational wave can traverse large amounts of dust without suffering a significant 
absorption by the medium. This suggested to us the idea that the long, filamentary sequences 
of galaxies seen in the large scale structure of the Universe may act as "optical" guides for the 
gravitational radiation that propagates in their interior: the gravitational wave from a source 
placed at one end of such a filament (or inside it) may get trapped by the large scale galactic 
gravitational field and consequently be able to traverse huge distances without undergoing 
the 1/R attenuation that characterizes propagation in vacuum, resulting in a possibly strong 
magnification of the source for an observer placed at the other end of the guide. The possibility 
of an analogous phenomenon for the electromagnetic radiation was studied in ref . flOtl . 

In this Section we shall study this phenomenon in an highly idealized situation. We shall 
consider a point source S of monochromatic gravitational waves, placed at a large distance 
from an observer O. Between them, there is a long cylinder C filled with dust. We assume 
that S, O and the cylinder are at rest and that S and O are roughly aligned with the axis 
of the cylinder. For simplicity, we imagine that the distribution of dust inside the cylinder is 
stationary, while the density of matter outside the cylinder is negligible. Moreover, we assume 
that the background gravitational field Qab, generated by the dust, is weak (but much stronger 
than the perturbation) and that the relative motion of the dust is much slower than the speed 
of light, which is generally the case in astrophysical situations. Under these conditions the 
background can be described using the Newtonian approximation and so there will be a quasi- 
Minkowskian coordinate system {x'^} = {ct, x, y, z] where the dust is practically at rest and 
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the gravitational field y^jy of tlie background can be expressed in terms of the Newtonian 
potential $ generated by the dust: 

2$\ / 2$\ 



9oo = -[l + -^j , gij = 1^1- 5ij , 

i,j = 1,2,3 (51) 

with 

2|^>| 

\9,.u\ < < 1 • (52) 

We assume that the time variation of the potential <I>, due to the motion of the dust, is so slow 
that it can be neglected and thus we take $ to be time- independent. We choose the origin and 
orientation of the Minkowskian coordinate system such that the axis of the cylinder C coincides 
with the z-axis. Its top and bottom faces have radius D and have z coordinates respectively 
equal to z = and z = L. The source S is placed outside the cylinder at the point Ps of 
coordinates {x = a, y = 0, z = —I), with / > and \a\ <^ D. 

According to the WKB approximation developed in the previous Section, the propagation 
of the wave emitted by S is described by the eikonal equation, Eq. (^). Since the background 
is time-independent, we can make for the eikonal S{x) the ansatz: 

S{t,x,y,z) = k {S{x,y,z) - ct} , (53) 

where k = 27ri/o/c, I'o being the frequency of the wave. Plugging this expression in the eikonal 
equation (^) and using for the background the expression in Eq. ( ^l]) we find 



^{ {d^Sf + {dySf + {d.Sf } + ^(x, y, z) = i . (54) 

While deriving Eq. (|5^), we have divided the eikonal equation by 2g^^ , expanding inverses and 
products of g^u to first order in $/c^. Equation ( |5^ ) has the same form of the eikonal equation 
describing a light- wave propagating in a medium with a refractive index n = 1 — 2^'/c'^. 

In order to study the propagation of the wave inside the cylinder, we shall not search for 
the solution of Eq. ( |5^ ) that describes a spherical wave originating from P5, but we shall 
rather focus our attention on the corresponding rays x{v). We observe that Eq. ([5^ ) has 
the same form as the Hamilton-Jacobi equation for a unit mass particle, moving in R3 with 
Hamiltonian: 

H{x\ k) = + + ll) + U{x, y, z) , (55) 

where 

U{x,y,z) = 2$/c2 . (56) 
Thus, the rays x{v) are solutions of the Hamilton's equations: 

dx^^dH_ dk^ _dH i = i 2 3 (57) 

dv dli ' dv dx'^ ' ' ' 

subject to the constraint H = 1/2. 
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We imagine that, before reaching the bottom face Sg of the cyhnder, the wave traverses 
a region of space where the gravitational field is negligible and so the incoming wave entering 
the cylinder will not differ appreciably from a spherical wave with center in Ps- If Q is a 
point of So of coordinates (Cj^jO), the eikonal 5(C,^, 0) will be roughly equal to the distance 
R = [(C — fl)^ + C'^ + P]^^^ of Q from the source. We shall focus our attentions on the rays that 
cross Sq in directions approximately parallel to the axis of the cylinder and this requires 

« 1. (58) 

The weak-field condition, ^/c'^{x,y,z) <C 1, ensures that the direction of these rays will con- 
tinue to be practically parallel to the z-axis all along the cylinder and this allows us to use the 
paraxial approximation to study their propagation. We start by taking the following Taylor 
expansion of the eikonal of the incoming wave: 

From it we compute the components IxiCO aiid ly{C,S,) of the gradient of S{x) on the points 
of Eq. The third component IziCO should be computed from the constraint H = 1/2 (of 
course, we have to pick up the positive solution for Iz, as our wave propagates in the positive 
z-direction) but due to the paraxial condition IxiCiO 1' ^yidO ^ 1; and to the weak-field 
condition <I>/c^ <^ 1, we can make the approximation IziC^C) = 1- In conclusion, we can say 
with a good approximation that on Sq the rays x{v) associated with the spherical wave from 
S are such, that 

x(0;C,6 = C, y(0;C,e) = e, ^(0;C,0 = o, 
UO; c, = iy{0; C, 6 = f ^.(0; C, = 1 , (60) 

where we have labelled each ray by the coordinates (C, the point, where it crosses So, and 
we have arbitrarily chosen for the value of the parameter v at that point. Thus, in order 
to see how the wave propagates inside the cylinder, in a small solid angle parallel to the axis, 
we just have to find the solutions of the Hamilton's equations ([57|), using Eqs. (^) as initial 
conditions. Consider now the case that the cylinder is much longer than thicker, L ^ D and 
that the distribution of dust inside it is independent of z so that, away from the caps, the 
potential $ will be only a function of {x,y). The z-components of Eqs. ( |57| ) are then trivial 
to integrate and give z{v) = v and lz{v) = 1. We can thus use z as a parameter in the 
place of V. As for the projections of Eqs. (^) in the (x,y) plane, they are formally identical 
to the equations of motion for a point particle on a plane, having a mass equal to one and 
moving in the potential U{x,y). It is obvious that if the potential U{x,y) has the shape of a 
two-dimensional well, there may exist a subset S* C So, such that for {(,, ^) G S*, the solutions 
{x{z; C)0) yi^'^ C)0) represent bounded motions and, when this happens, it means that the 
rays that cross S* will get trapped inside the cylinder. We thus have that the gravitational 
field of the cylinder may effectively act as a sort of waveguide for the gravitational radiation. 
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In order to examine this phenomenon in greater detail, let us make the simplifying assumption 
that the density pQ of the dust is uniform. Then ?7(x, y) inside the cylinder is equal to the 
potential of a harmonic oscillator: 



U{x,y) = \oj'^{x^ +y'^) for (x^ + y^) < ^ 



where 



UJ 



(61) 



(62) 



If D is infinite, the motions (x(z; CiC)-, vi^'i CO) are all bounded and thus all the rays will 
be trapped by the cylinder. Upon solving the Hamilton's equations, we get 



1 



cos{uJz) H sin(ujz] 



uil 



cos(wz) H sin(a;2;) 



x{z-. CO = c 
y{z; CO = e 

Let now Lfoc be the solution of the transcendental equation: 

tan(LL'Lfoc) = —luj it/2 < wLfoc < vr . 



sin(a;z). 



(63) 



(64) 



Assume that L > Lfoc and let be the integer such that NLfoc < L < {N + l)Lfoc- We see 
from Eqs. (p3|) that the rays are focussed at the focal points -F„ of coordinates 



Pn = {-(-l)"acos(wLfoc 

n = l,2,---iV 



vr , 



0, Lfoc+ (n- 1)-} , 

UJ 



(65) 



This means that inside the cylinder there will form a series of alternatively inverted and right 
images of the source. We shall now compute the specific luminosity of the image at F^, 
assuming that the source S has a specific luminosity £5. According to the definition of Cs, 
and assuming that S radiates isotropically, the number dNs of gravitons with energy in the 
range hdiVs emitted by S during the proper time interval drs in the solid angle dils around 
the z-axis is equal to 



dNs = drs dils dtos 



Attljs 



(66) 



These gravitons form a narrow bundle. According to what we said above, they will be focussed 
at Fn, and we let di^n be the solid angle they span at F„. By the graviton number conservation, 
Eq. (|50|), dNs must be equal to the number of gravitons dNn emitted by the image during the 
corresponding time interval dr^, in the interval of frequency dun in the solid angle dO^: 



dNs = dNn = dTn dQn duJn 



(67) 



Since the source and the image are at rest, drs = dTn-, ujs = and dus = dwn- Thus, we 
have 



d^c 

dClr 



(68) 
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In view Eqs. (p3|), it is easy to verify that 



cos^(cuLfoJ(l + u;^/^)^ . (69) 
alls 

If <C 1 we see from Eq. ( |64D that LfQ^ ~ tt/uj — I. Equation ( |69D then imphes Hn — 
Cs{l —uj'^P) ~ Cs, which shows that the images have practically the same specific luminosity 
of the source. Prom the point of view of the observer O, it will be as if the "waveguide" had 
driven the source from Ps to F/v- Being much closer to him than the real source, the image 
at -F/v will obviously appear to him much brighter. 

So far, we have proceeded as if the radius of the cylinder had been infinite. The effect 
of D being finite is that in general not all the rays entering the cylinder will be focussed at 
F/v. Outside the cylinder, the strength of the gravitational field decreases rapidly as one goes 
farther from its axis and so in general a ray that escapes out of it will not be focussed at Fjy. 
This suggests that a simple criterion for a ray to be focussed is that it should always remain 
inside the cylinder while traversing it: + < for < z < L. Assuming again that 

Lol <^ 1, we see from Eqs. (^) that this condition is satisfied only by the gravitons emitted by 
5 in a solid angle A^ls around the direction of the z-axis of order: 

Ans « -^vr . (70) 

Since, for col <^ 1, this is also roughly equal to the angle An^r spanned by these rays at the 
focal point F/v, Ail^ provides a measure of the degree of alignment of O with the axis of the 
cylinder needed for O to observe the image at F/v. 

In order to have an estimate of the above numbers in realistic situations, let us imagine 
that our waveguide is constituted by a sequence of galaxies forming a filamentary structure of 
the kind observed in the large scale structure of the Universe. Taking for D the typical size of 
a galaxy, D ^ 10^ pc, and for po the typical density of matter inside a galaxy, ~ 5 • 10~^^ 
g/cm^, we find that the weak-field approximation for the background, the second inequality in 
Eq. (^), is satisfied, because near the edge of the cylinder, where the galactic gravitational 
field is strongest, 2<I>/c^ ~ a;^F>^/2 ~ 10~^. Plugging the above value for po in Eqs. ( |6T|) 
and (|6^), we get that Lfoc ~ 10^ Mpc, which is the typical scale of the large structures of 
the Universe and represents also the expected distance of coalescence events, as was discussed 
in the Introduction. The corresponding angular opening AQs of the bundle of rays that are 
efficiently trapped is of the order of 1'. 

In the above discussion, we have treated our waveguide as if the distribution of matter 
inside it was homogeneous. Obviously this is not what happens inside the galaxies, where a 
significant fraction of matter is concentrated in compact bodies like the stars. |^ The large scale 

^Considering the recent results of microlensing collaborations like MACHO, OGLE, DUO, and EROS, it is 
very likely that a large amount of the matter inside galaxies is built up of baryonic compact objects weakly 
interacting with electromagnetic radiation like planets, brown dwarfs or black holes so that all the following 
considerations apply also to these galactic components [p^ , po[ , pl|] , . 
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galactic field is distorted on small scales by the local field of nearby stars and we have to analyze 
the effect of these local disturbances on the focussing properties of the waveguide. A simple 
order-of-magnitude estimate shows that in fact it should be negligible. We start by computing 
the distance d{x) such that a ray, while traversing the cylinder, passes with probability a; at a 
distance smaller than d{x) from a star. For this purpose, consider the tube T of radius d{x) 
around the ray. Assuming that the ray is roughly parallel to the axis of the cylinder and that 
the stars are distributed in the cylinder with uniform density n^,, the probability of finding a 
star inside T is 

X = 1 - (1 - d2(^)/£)2)n.^D2L ^ n,7rd^{x)L . (71) 

This formula is easy to understand, since {l — d'^{x)/D^)"'*'^^^^ is the probability that throwing 
at random n^firD^L stars inside the cylinder no one of them falls inside T. Now, since the 
deflection a{x) suffered by a ray passing at a distance d{x) from a star with mass M is 
approximately a{x) ~ AGM / (?d{x), we have 



Taking L ~ 10^ Mpc, Mn* ~ 0.12 Mqpc^^ p^ ] and assuming that the stars have a mass of 
the order of the mass of the sun Mq , we get 

a{x) ~ 2 X 10"^^ arcsec. (73) 



For X = 0.01, the deflection angle is of order of 10~^ arcsec. This means that, while traversing 
the cylinder from one end to the other, a ray has a probability of 1% of being deflected by a 
star by an angle of approximately one thousandth of arcsecond. Now consider the bundle of 
rays emitted by 5" in a solid angle of opening AJ^^. While traversing the cylinder, the cross 
section of this bundle will expand up to a point where it becomes comparable to the cross 
section of the cylinder and then, after traversing a distance of order Lfoc, it will shrink to a 
point. The rays in the bundle will thus suffer an average deflection of order D /L^oc — I'l which 
is much greater than the deflection due to the occasional passages near a star as estimated 
above. We thus conclude that the stars have a negligible effect on the propagation of the rays 
inside the cylinder. 



4 Discussion and conclusions 

In this paper, we have shown that a filamentary distribution of dust may act as a waveguide 
for gravitational radiation. We derived this result studying the short wave-length limit of 
the linearized Einstein's equations that describe the propagation of a weak gravitational wave 
through dust. In this approximation, a gravitational wave behaves in essentially the same 
way as an electromagnetic wave propagating in vacuum: in particular, the gravitational wave 
propagates along geodesies of the background gravitational field generated by the dust. The 
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dust itself behaves as a non-dispersive and non-dissipative medium, preserving the polarization 
state of the gravitational wave. In this geometrical optic limit, the waveguiding effect has a very 
simple explanation: the gravitational field generated by a filament of dust traps in its interior 
a paraxial bundle of rays, with the result that, while traversing the filament, the amplitude h 
of the gravitational wave does not undergo the 1 /R attenuation characterizing the propagation 
in vacuum of spherical waves. We showed also that for a roughly uniform density of dust, the 
guide produces real images of the source, having approximately the same absolute "luminosity" 
as the real source. In some sense, the waveguide "draws" the source closer to the observer: if 
the true distance to the source is R, its image brightness will correspond to that of a similar 
source at the closer distance 

Reff = R-l-L, (74) 

where L is the length of the filament and / is the distance of the source from the top of the 
filament (see also [|TC|]). 

Even though the effect described in this paper is basically a gravitational lensing effect, it 
is worth pointing out that it differs in some important points from the standard gravitational 
lensing situations. When discussing the gravitational lensing of light, it is normally sufficient 
to treat the deflector as a thin lens; moreover it is usually assumed that the gravitational field 
of the deflector can be treated using the newtonian limit, which requires the impact parameter 
ro to be much larger than the gravitational rg radius of the lens vq ^ rg. This is not a very 
stringent constraint, because of the smallness of r^; a much more stringent condition is that 
the light ray should pass close enough to the lens, for the deflection to be sensible, but not so 
close to cross it, because this would easily lead to the absorption or scattering of the light ray 
by the matter constituting the lens. In order to obtain the waveguiding effect described in this 
paper, we need to reverse two of the three above conditions: while we still assume the weak 
field condition, which is well respected in realistic astrophysical situations, an efficient trapping 
of the gravitational wave requires the lens to be very thick (in fact our filamentary lenses are 
much thicker than wider), and that the wave passes well inside the distribution of matter, in the 
region where the newtonian potential well is deeper and thus more confining. These conditions 
do not represent a problem since gravitational waves, differently from electromagnetic waves 
p4[,p^], travel nearly unscathed through all forms and amounts of intervening matter 
and can thus traverse our long filaments without being significantly absorbed or scattered. 
Finally, we should mention that even though in this paper we have limited the discussion of the 
waveguiding effect to a dust distribution with the shape of a filament, more general geometries, 
planar for example, analogous to those used for optical waveguides, could be considered. 

If the waveguiding effect described in this paper could be really observed, it might signif- 
icantly lower the sensitivity threshold required to successfully detect a gravitational wave. In 
order to address this issue several questions have to be answered. The first one is: what is 
the probability that waveguides for gravitational radiation exist? Furthermore, what is the 
probability that a significant fraction of the flux from sources of gravitational radiation will 
get trapped? 
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The first point seems, from our point of view, very likely. The observations in the last 
twenty years tell us that filamentary or slab-like structures are quite common at large scales 
(see, e.g. ||2^,[^,|^,Q). Furthermore, their sizes (~ 10^ 10^ Mpc) are comparable 
with the focalization length of waveguides, as we have shown above. The hypotheses of their 
origin, starting from some fundamental theory and in connection with primordial structures as 
cosmic strings, domain walls and textures are among the most investigated topics of theoretical 
cosmology pi , pi] , P , H . 



The second point depends on the efficiency of production of gravitational radiation. As we 
said in the Introduction, for rates of the order of a few events per year, the expected distance 
of high frequency sources like coalescing neutron-stars, should be in the range of 10^ 10^ 
Mpc, which is also the order of magnitude of the focalization length of our waveguides. If the 
probability that a wave from such a source gets trapped by a waveguide is large enough, this 
effect might be seen already by the VIRGO and LIGO interferometers in the next few years. 
Since the filamentary structures that constitute our waveguides extend over cosmological dis- 
tances, the waveguiding effect may have an important impact on the detection of gravitational 
waves in the low-frequency band that will be the target of the LISA interferometer and the 
candidates for such a kind of sources are quite numerous. First of all we have to stress the fact 
that, obviously, all the sources we are considering have extragalactic origin and can be sited 
in distant galaxies. This fact implies that our sources may have high redshifts so that in the 
frequency and in all the quantities characterizing a gravitational wave we have to include a 
factor (l + z). A large amount of candidates sources could be formed by binary evolved stars as 
double systems of white dwarfs, black holes and white dwarfs, or black holes and black holes. 
Such structures are very common in the Universe, they are present in all galaxies and, if they 
are tight bound, the gravitational radiation emission fluxes are conspicuous. Another class of 
objects could be quasars if in their core a very massive black hole is present and massive 
objects fall into it. Finally, waves from primordial phase transitions and from cosmic strings 
p2[ could be good candidate sources. 
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